International  Journal  of  Power  Electronics  and  Drive  System  (IJPEDS) 

Vol.  10,  No.  1,  March  2019,  pp.  137-150 

ISSN:  2088-8694,001:  10.1 1591/ijpeds.vl0.il.ppl37-150 


□  137 


Output  feedback  nonlinear  control  of  three-phase 
grid-connected  PV  generator 


A.  Yahya,  H.  El  Fadil,  M.  Oulcaid 

Systems  Engineering  Laboratory  (LGS),  National  Sehool  of  Applied  Seienees  (ENSA),  Ibn  Tofail  University,  Moroeeo 


Article  Info 


ABSTRACT 


Article  history: 

Received  May  16,  2018 
Revised  Aug  18,  2018 
Accepted  Sep  3,  2018 


Keywords: 

MPPT 

Nonlinear  control 
Nonlinear  observer 
Output  feedback  control 
Renewable  energy 
Stability  analysis 


This  paper  addresses  the  problem  of  eontrolling  three-phase  grid  eonneeted 
PV  system  involving  a  PV  arrays,  a  voltage  souree  inverter,  a  grid  filter  and 
an  eleetrie  grid.  This  paper  presents  three  main  eontrol  objeetives:  i)  ensuring 
the  Maximum  power  point  traeking  (MPPT)  in  the  side  of  PV  panels,  ii) 
guaranteeing  a  power  faetor  unit  in  the  side  of  the  grid,  iii)  and  ensuring  the 
asymptotie  stability  of  the  elosed  loop  system.  Interestingly,  the  present 
study  features  the  aehievement  of  the  above  energetie  goal  without  resorting 
to  sensors  of  eurrents  of  the  grid.  To  this  end,  an  output- feedbaek  eontrol 
strategy  eombining  a  state  observer  and  a  nonlinear  eontrol  laws  is 
developed.  The  proposed  output-feedbaek  eontrol  strategy  is  baeked  by  a 
formal  analysis  showing  that  all  eontrol  objeetives  are  aetually  aehieved. 
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1.  INTRODUCTION 

It  was  very  clear  from  recent  studies  and  documentation  that  fossil  fuels  would  last  only  a  few  more 
decades.  The  cost  of  fossil  fuels  has  become  a  major  challenge  for  all  of  humanity.  Not  only  the  economic 
value  but  also  the  environmental  impacts  of  fossil  fuels  have  clearly  pushed  us  towards  alternatives.  The 
biggest  alternatives  that  can  really  make  a  difference  to  sustainability,  such  as  reducing  greenhouse  gas 
emissions  and  the  long-term  economy,  are  renewable  energy  sources  (RES)  such  as  wind  and  solar  energy 
[1],  [2].  PV  system  is  increasing  as  a  renewable  source  due  to  its  advantages  of  little  maintenance,  absence  of 
moving  mechanical  parts,  no  noise  and  no  pollutant  emission  [3].  Furthermore,  its  cost  is  decreasing  over  the 
next  ten  years,  while  the  deployment  of  PV  systems  continues  to  increase  rapidly.  With  increasing  PV 
penetration  on  the  grid,  eventually  reaching  hundreds  of  gigawatts  (GW)  of  interconnected  capacity,  a  variety 
of  methods  must  be  considered  and  implemented  at  different  scales  for  a  reliable  and  cost-effective 
connection  into  power  grid. 

Different  control  strategy  for  three-phase  grid  connected  of  PV  modules  has  been  largely  dealt  with 
in  the  specialist  literature  in  the  last  few  years  [4]-[9].  Nevertheless,  good  integration  of  medium  or  large  PV 
system  in  the  grid  may  therefore  require  additional  functionality  from  the  inverter,  such  as  control  of  reactive 
power.  Moreover,  the  increase  in  the  average  size  of  a  PV  system  may  lead  to  new  strategies  such  as 
eliminating  the  DC-DC  converter,  which  is  usually  placed  between  the  PV  generator  and  inverter,  and 
moving  the  MPPT  to  the  inverter,  which  leads  to  increased  simplicity,  overall  efficiency  and  a  cost  reduction. 
These  two  features  are  present  in  the  three-phase  inverter  that  is  presented  in  this  paper,  with  the  addition  of  a 
Perturb  and  Observe  (P&O)  MPPT  algorithm. 
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The  present  paper  is  then  foeusing  on  the  problem  of  eontrolling  three  phase  grid-eonneeted  PV 
power  generation  systems.  The  eontrol  objeetives  are  threefold:  (i)  global  asymptotie  stability  of  the  whole 
elosed-loop  eontrol  system;  (ii)  aehievement  of  the  MPPT  for  the  PV  array;  and  (iii)  ensuring  a  grid 
eonneetion  with  unity  Power  Faetor  (PF).  These  objeetives  should  be  aehieved  in  spite  the  elimatie  variables 
(temperature  and  radiation)  ehanges.  To  this  end,  a  nonlinear  eontroller  is  developed  using  Lyapunov  design 
teehnique.  A  theoretieal  analysis  is  developed  to  show  that  the  eontroller  aetually  meets  its  objeetives  a  faet 
that  is  eonfirmed  by  simulation. 

The  paper  is  organized  as  follows:  the  three-phase  grid  eonneeted  PV  system  is  deseribed  and 
modeled  in  Seetion  2.  Seetion  3  is  devoted  to  eontroller  design  and  analysis.  The  eontroller  traeking 
performanees  are  illustrated  by  numerieal  simulation  in  Seetion  4. 


2.  RESEARCH  METHOD 

2.1.  System  overview 

The  main  eireuit  of  three-phase  grid-eonneeted  photo voltaie  system  as  shown  in  Figure  1 .  It  eonsists 
of  a  PV  arrays;  a  DC  link  eapaeitor  C;  a  three  phase  inverter  (ineluding  six  power  semieonduetors)  that  is 
based  upon  to  ensure  a  DC-AC  power  eonversion  and  unity  power  faetor;  an  induetor  filter  L  with  a  series 
resistanee  r,  and  an  eleetrie  grid.  The  eontrol  inputs  of  the  system  are  a  PWM  signals  u^and  taking 
values  in  the  set  {0,1}.  The  grid  voltages  Vga,  Vgjj  and  Vg^  eonstitute  a  three-phase  balaneed  system. 

2.2.  PV  Array  model 

An  equivalent  eireuit  for  a  PV  eell  as  shown  in  Figure  2.  Its  eurrent  eharaeteristie  ean  be  found  in 
many  plaees  [10]-[13]  and  presents  the  following  expression. 


PV  Generator  DC  Link  Inverter 

Figure  1.  Three-phase  grid  eonneeted  PV  system 
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The  meaning  and  typieal  values  of  the  parameters  given  by  (1)  and  (2)  ean  be  found  in  many  plaees 
(see  e.g.  [14],  [15],  [16]).  A  is  diode  ideal  faetor,  k  is  Boltzmann  eonstant ,  T  is  temperature  on  absolute  seale 
in  K,  q  is  eleetron  eharge  and  k  is  the  radiation  in  kW/m^,  Iphr  is  the  short-eireuit  eurrent  at  298  K  and  1 
kW/m^,  Ki  =  0.0017A/K  is  the  eurrent  temperature  eoeffieient  at  Iphr,  Ego  is  the  band  gap  for  silieon,  is 
referenee  temperature,  Latr  is  eell  saturation  eurrent  at  Tr.  PV  array  eonsists  of  Ns  eells  in  series  formed  the 
panel  and  of  Np  panels  in  parallel  aeeording  to  the  rated  power  required.  The  output  voltage  and  eurrent  ean 
be  given  by  the  following  (3)  and  (4). 
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Figure  2.  Solar  cell  circuit  diagram 


The  photovoltaic  generator  considered  in  this  paper  consists  of  several  NU-183E1  modules.  The 
corresponding  electrical  characteristics  of  PV  modules  are  shown  in  Table  1.  The  associated  power- voltage 
(P-V)  characteristics  under  changing  climatic  conditions  (temperature  and  radiation)  are  shown  in  Figures  3 
and  4.  This  highlights  the  Maximum  Power  Point  (MPP)  Ml  to  M5,  whose  coordinates  are  given  in  Table  3. 
The  data  in  Table  1  to  Table  3  will  be  used  for  simulation. 

Table  2  shows  the  main  characteristics  of  the  PV  array,  designed  using  Sharp  NU-183E1  modules 
connected  in  a  proper  series-parallel,  making  up  a  peak  installed  power  of  71  kW.  As  there  is  no  DC-DC 
converter  between  the  PV  generator  and  the  inverter,  the  PV  array  configuration  should  be  chosen  such  that 
the  output  voltage  of  the  photovoltaic  generator  is  adapted  to  the  requirements  of  the  inverter.  In  this  case  a 
380V  grid  has  been  chosen,  so  the  inverter  would  need  at  least  570V  DC  bus  in  order  to  be  able  to  operate 
correctly.  The  minimum  number  of  modules  connected  in  series  should  be  determined  by  the  value  of  the 
minimum  DC  bus  voltage  and  the  worst-case  climatic  conditions.  The  PV  array  was  found  to  require  28 
series  connected  modules  per  string. 


Table  1.  Electrical  specifications  for  the  solar  module  NU-183E1 


Parameter 

Symbol 

Value 

Maximum  Power 

Pm 

183W 

Short  eireuit  eurrent 

Iscr 

8.48A 

Open  eireuit  voltage 

Voc 

30.1V 

Maximum  power  voltage 

Vm 

23.9V 

Maximum  power  current 

Im 

7.66A 

Number  of  parallel  modules 

Np 

1 

Number  of  series  modules 

Ns 

48 

Figure  3.  (P-V)  characteristics  of  The  PV  Generator 
(NP=14  And  NS=28)  with  constant  temperature  and 
varying  radiation 


Figure  4.  (P-V)  characteristics  of  The  PV  Generator 
(NP=14  And  NS=28)  with  constant  radiation  and 
varying  temperature 
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Table  2  :  PV  Array  specifications  using  sharp  NU- 


183E1 


Parameter 

Symbol 

Value 

Total  peak  power 

Ptm 

71  kW 

Number  of  series  strings 

NS 

28 

Number  of  parallel 

NP 

14 

Number  of  PV  panels 

N 

392 

Voltage  in  maximum  power 

Vm 

664V 

Current  peak 

Im 

107A 

Table  3.  Maximum  power  points  (MPP):  in  Figure  3 

_ and  Figure  4 _ 

MPP _ Vm[  V] _ Pm  [  KW  ] 


Ml 

664.2 

71 

M2 

661 

57.2 

M3 

648.8 

35.1 

M4 

664.2 

71.5 

M5 

584 

61.7 

2.3.  Modeling  of  three-phase  grid-connected  PV  system 

The  state-space  model  of  a  three-phase  grid-connected  photovoltaic  system  shown  in  Figure  1  can 
be  obtained  by  the  dynamic  equations  described  in  (5a)-(5d): 
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Applying  the  Concordia  transformation  to  (5a-d),  the  instantaneous  model  in  stationary  coordinates  is  givens 
as  (6a)-(6b): 
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And  the  transformation  matrix  is  given  by: 
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(6b) 

(6c) 

(7a) 


(7b) 
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According  to  the  three-phase  instantaneous  active  and  reactive  power  theory,  the  instantaneous 
active  power  P  and  the  instantaneous  reactive  power  Q  can  be  derived  in  three-phase  abc  coordinates  for  the 
three-phase  photovoltaic  grid-connected  inverter: 


P  =  VgJa+Vgbh+'^Jc 

(8a) 

Q  =  [(v  -  V  X  +  (v  -  V  X  +  (v  -  V  X 

(8b) 

In  the  af  coordinates,  the  instantaneous  active  power  P  and  the  instantaneous  reactive  power  Q  can 
be  expressed  as  following: 

P  = 

(9a) 

e= 

(9b) 

3.  OUTPUT  FEEDBACK  CONTROLLER  DESIGN 

The  model  (6a-c)  is  useful  to  build  an  accurate  simulator  for  the  studied  system  however  it  is  not 
adequate  to  elaborate  a  continuous  controller  as  it  involves  a  binary  input  u^  and  Up  .  For  control  design 

purpose,  it  is  more  convenient  to  consider  the  following  averaged  model,  obtained  by  averaging  the  model 
(6a-c)  over  one  switching  period  [17]. 

r  X.  1 

(10a) 

r  X.  1 

(10b) 

X}  =  +  Hl)X^+  V 

(10c) 

where  ,  X2  ,  X3  ,  ju^  and  ju^  denote  the  average  values  of,  respectively  and 

Up  .  The  signals  ju^  and  jUp  ,  called  duty  ratios  which  belong  to  [0,l] ,  are  considered  as  the  inputs  of  the 
system  (lOa-c). 

On  the  basis  of  (lOa-c),  the  next  two  subsections  will  be  devoted  to  the  observer  design  and  the 
nonlinear  controller  design. 

3.1.  Observer  design 

Let  us  now  consider  in  (lOa-c)  that  is  the  only  measurable  variable.  Then  the  following  nonlinear 
observer  is  proposed: 


A  r  ^  X.  1 

(11a) 

;  r  ,  X,  1 

X2=--x^+—Hp--y^p 

(11b) 

X3=-fMji  +MfiA+fp^  +  ^(x3 -X3) 

(11c) 

Where  ,  X2 ,  X3  represent  the  estimates  of  the  state  variables  and  T  >  0  being  the  observer  design 
parameter. 

Let  us  introduce  the  estimation  errors  =  x^  -x^ ,  Z2  =  X2  -X2  and  Z3  =  X3  -X3 .  Then,  from  (lOa-c) 
and  (lla-c),  one  has: 
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V  7 

(12a) 

T  Z 

4=--Z2+Y/^/, 

(12b) 

=-^(//„Zi+//^Z2)-/lZ3 

(12c) 

From  the  error  system  (12a-e),  one  ean  state  the  following: 

Proposition:  Consider  the  estimation  error  system  (12a-e),  obtained  by  eombining  the  system  (lOa-e)  and  the 
nonlinear  observer  (1  la-e).  Then,  3  Tq  ,  V  T  >  Tq  ,  whatever  the  initial  eonditions,  the  state  estimation  error 

z  =  [zp  Z2,  Z3  ]  ^  eon  verges  exponentially  to  zero. 

Proof:  Consider  the  following  quadratie  Lyapunov  funetion: 

(13) 

Its  derivative,  using  (12a-e),  ean  be  obtained  as  follows: 


With 


-  Xzl  -1-  -1-  /^2^2^3 


1  1 


1  1 


L  C 


(14a) 


(14b) 


Using  Young’s  Inequality  (14a)  beeomes 


lA 
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2s^ 


2 


+  A 


1 


(15) 


where  and  ^2  being  any  real  positive  eonstants. 

As  ju^  G  [0,1]  and  e  [0,l]  (duty  ratio  funetions),  it  follows  that  \p^  <  P^  and  \p^  <  P^  where 


Po  = 


1 

L 


1 

C 


From  (15)  beeomes,  using  (16) 
<  —k^Z^  “^2^2  ~^3^3 


(16) 


(17a) 


Where 


K='t-P-  k,=^-^{s^  +  s^) 


L  2s^ 


L  2s. 


(17b) 


Let  s.  and  S.  to  be  ehosen  as  follows:  s. 


and 

2r 


£.  > 


.  Then,  the  observer  design  parameter  ^ 


2r 


ean  be  ehosen  so  that  X>  where: 
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\=^{s,+s,)  (18) 

Therefore  >  0 ,  A:2  >  0  and  A:3  >  0  .  It  follows  that  (17a)  ean  be  rewritten  as  follows: 

(19) 

where 

k  =  2  min(A:j  ,k2,k^)  (20) 

From  (19),  it  is  easy  to  show  that: 

(21) 


whieh  yields  that  V^(t)  is  exponentially  vanishing.  It  follows  that  limz(^)  =  0  ,  whieh  in  turn  shows  that  the 

>00 

estimates  eon  verge  toward  their  true  values.  This  ends  the  proof  of  the  Proposition.  In  the  next  subseetion,  we 
foeus  on  elaborating  a  eontroller  that  stabilizes  the  system. 

3.2.  Nonlinear  controller  design 

With  the  aim  of  design  an  appropriate  eontrol  for  the  model  (lOa-e)  deseribed  in  previous  seetion, 
the  eontrol  objeetives,  the  eontrol  design  and  stability  analysis  will  be  investigated  in  this  Seetion,  taking  into 
aeeount  the  nonlinear  feature  and  the  multi-input  multi-output  (MIMO)  aspeet  of  the  system.  In  order  to 
define  the  eontrol  strategy,  the  first  step  is  to  establish  eontrol  objeetives,  whieh  ean  be  summarized  as 
follows: 

a.  Maximum  power  point  traeking  (MPPT)  of  PV  arrays, 

b.  Unity  power  faetor  (PF)  in  the  grid, 

e.  Asymptotie  stability  of  the  whole  system. 

It  is  worth  noting  that  all  objeetives  must  be  aehieved  without  sensing  all  variables.  The  PV  voltage 
X3  is  the  only  measurable  variable. 

The  first  eontrol  objeetive  is  to  enforee  the  real  power  P  to  traek  the  maximum  power  point  Pm>  It’s 
already  point  out;  that  this  power  ean  be  eontrolled  by  the  a-axis  eurrent  4  and  fi-axis  eurrent  ip  (see  (9a)).  In 
this  paper  the  MPPT  algorithm  based  on  the  Perturb  and  Observe  (P&O)  teehnique  [18]  is  resorted  to 

generate  the  eoeffieient  ^  (involved  in  (22))  so  that  the  aetive  power  P  traeks  its  maximum  value  i.e. 

P  =  Pm  • 

The  seeond  eontrol  objeetive  means  that  the  grid  eurrents,  i^ ,  i^,  and  i^  should  be  sinusoidal  and  in 
phase  with  the  AC  grid  voltage  and  respeetively.  To  this  end  the  reaetive  power  have  to  be  null. 

To  aehieve  this  objeetives  it  suffiees  to  enforee  the  a-axis  eurrent  i^  and  p-axis  eurrent  i^  to  traek 
referenee  signals,  say  x^  and  xj ,  of  the  following  forms: 


xl=0-v^^ 

(22) 

x\=0-v^l, 

(23) 

With  ^  is  any  real  positive  parameter  (although  transient  time-variations  are  allowed).  Its  generation  will  be 
seen  later  using  the  MPPT  algorithm  (see  seetion  4).  Onee  the  eontrol  objeetives  are  elearly  defined,  as  the 
MIMO  system  is  highly  nonlinear,  a  Lyapunov  based  nonlinear  eontrol  is  proposed  [17].  Given  the  observer 

guarantees  that  the  errors  ^  1  1  and  222  converge  to  zero,  the  following  eontroller  design  will 

be  based  on  the  estimate  variables  instead  of  their  true  values.  Therefore,  the  following  errors  are  introdueed: 

z^=x^-xl  (24) 
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/V  * 

Z5  —  2C2  2C2 

(25) 

In  order  to  aehieve  the  objeetives:  MPPT  and  power  faetor  unit,  one  ean  seek  that  the  errors  Z4  and  z^  are 
vanishing.  To  these  ends  the  dynamies  of  Z4  and  z^have  to  be  elearly  defined.  Deriving  (24)  and  (25),  it 
follows  from  (11a)  and  (11b)  that: 

^4  =  “7^1  +  “  7  V 

(26) 

r  ^  X.  1  .* 

^5=--^2+YP^--V-^2 

(27) 

The  goal,  now,  is  to  make  Z4  and  z^  exponentially  vanishing  by  enforeing  its  derivatives 
behave  as  follows: 

Z4  and  Zj  to 

Z4  =  -A  sgn(zj  -  z^{T)dT 

(28) 

=  -Pi  sgn(z5 )  -  ^2  J  Z5  i'^)d  z- 

(29) 

where  >  0  ,  >02  >  0  ,  >  0  and  4  >  b  being  design  parameters  and  is  a  signum  funetion.  It  is  worth  noting 
that  the  integral  aetions  are  introdueed  in  (28)  and  (29)  to  allow  a  good  robustness  of  the  eontroller  against 
unmodeled  dynamies  and  perturbations.  Combining  (26)  and  (28)  the  first  eontrol  law  ju^  is  obtained: 

Pa  ASgn(Z4)  +  -2x,  +jv^^+x’ z^{T)dT^ 

(30) 

Finally,  combining  (27)  and  (29),  the  second  control  law  /Up  is  also  obtained 

Pii=y{^P^  Sgn(z5 )  +  ^ Xj  +  +  xj  -  #2  j Z5 (r)rfr  j 

(31) 

Since  the  eontrol  laws  are  elearly  defined,  the  eoneem  now  is  to  investigate  the  eonvergenee  of  the  errors  Z4 
and  ^5 .  To  this  end  the  following  quadratie  Lyapunov  funetion  is  eonsidered 


Vr  +-^[|z4(r)^/r|  + -y  [|z5(r)^/r]^ 


1  2,1, 

2^"^  2^"  2' 


Its  derivative  is  obtained  as  follows  Its  derivative  is  obtained  as  follows: 

Vc  =24^4+  Z5Z5  +  ^1^4  [j  Z4  r  +  4^5  [j  Z5 

whieh,  using  (28)  and  (29),  gives: 

“  “Pi  1^4 1  “  P2 1^5 1 


(32) 


(33) 


(34) 


As  Vc  is  positive  definite  funetion  and  its  derivative  (34)  is  negative  definite  it  follows  that  the  equilibrium 

(z4,Z5)  =  (0,0)  is  globally  asymptotieally  stable  [18].  Whieh,  in  turn,  gives  lim(z4 (^),  z^ (^))  ^  (0,0) 

>00 

The  main  result  of  the  proposed  output  feedback  controller  is  summarized  in  the  following  theorem. 
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Theorem:  Consider  the  closed-loop  system  consisting  of  the  controlled  system  of  Figure  1  represented  by  its 
nonlinear  model  (lOa-c),  the  nonlinear  observer  (lla-c)  and  the  controller  composed  of  the  control  laws  (30) 
and  (31).  Then,  one  has: 

a.  The  closed  loop  system  is  GAS.  It  follows  that  all  closed  loop  signals  are  bounded. 

b.  The  estimation  errors  z  =  [z^ ,  Z2 ,  Z3  ]  ^  converges  exponentially  to  zero. 

c.  The  tracking  errors  and  converge  to  zero  implying  MPPT  achievement  and  power  factor  unit. 


Proof:  let  us  consider  the  following  quadratic  Lyapunov  function: 


V  =  V  +V  = 

^  2 


^  +  ^^Z^(T)dT 


If^ 


+  ^[h5(z'Vz' 


(35) 


Its  derivate,  using  (19)  and  (34),  can  be  obtained  as  follows: 
V  =  Vo  +  Vc<  -|(zf  +  +  zl)-  p,\z,\- p,\z,\ 


(36) 


Which  clearly  shows  that,  the  equilibrium  z  =  0  of  the  closed  loop  system  with  the  state  error  vector 
z  =  [zj,Z2,Z3,Z4,Z5]^  is  globally  asymptotically  stable.  It  follows  that  all  errors  are  vanishing.  This  ends  the 
proof  of  theorem.  The  next  section  devoted  to  the  performances  evaluation  of  the  proposed  output 
feedback  controller. 


4.  SIMULATION  RESULTS 

The  theoretical  performances  described  by  the  theorem  of  an  output  feedback  controller,  including 
the  control  laws  (30-31)  and  the  nonlinear  observer  (lla-c),  designed  in  Section  3,  are  now  illustrated  by 
simulation.  The  experimental  setup,  described  by  Figure  5,  is  simulated  using  MATLAB/SIMULINK.  The 
characteristics  of  the  controlled  system  are  listed  in  Table  4.  Note  that  the  controlled  system  is 
simulated  using. 


MPPT 

l^rturb 

Observe 


IK7AC 

InvfiOer 


TmTT 

^flJ 

f  1  i 


apiabt 


i.t 


Conirdicr 

liwt(ja)  Mfid 

t 


Nanlinesr 


t  t 

Figure  5.  Simulation  bench  of  the  proposed  three-phase  grid  connected  system 


The  instantaneous  three  phase  model  given  by  (5a-d).The  model  (lOa-c)  in  a-p  axis  is  only  used  in 
the  controller  design.  The  design  parameters  of  the  controller  are  given  values  of  Table  5.  These  parameters 
have  been  selected  using  a  Trial-and-error’  search  method  and  proved  to  be  suitable.  It  worth  noting  that  the 
parameter  0  involved  in  (22)  and  (23)  is  generated  using  Perturb  and  Observe  algorithm  with  the  block- 
diagram  illustrated  by  Figure  6.  The  resulting  closed  loop  control  performances  are  illustrated  by  Figure  7 
to  Figure  13. 
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Table  4.  Characteristics  of  controlled  system 


Parameter 

Symbol 

Value 

PV  array 

PV  power 

40  kW 

DC  link  eapaeitor 

C 

3300pF 

Grid  filter  induetor 

L 

r 

3mH 

0.2Q 

PWM 

Switehing  frequeney 

lOkHz 

Grid 

AC  souree 

Line  frequency 

220V 

50Hz 

(Td) 


Figure  6.  P&O  algorithm  implementation  in  Matlab/Simulink  software 


Table  5.  Controller  parameters 


Parameter 

Symbol 

Value 

A 

10^ 

Pi 

4x10^ 

Design  parameters 

0.1 

0.1 

/I 

5x10^ 

Delay  time  Td 

10"'* 

P&O  algorithm  parameters 

Step  value  k 

U.3 

4.1.  Radiation  change  effect 

Figure  7  shows  the  perfect  MPPT  in  the  presence  of  radiation  step  changes  meanwhile,  the 
temperature  is  kept  constant,  equal  to  298.15K(25°C).  The  simulated  radiation  profile  is  as  follow:  a  first  step 
change  is  performed  between  500  and  lOOOW/m^  at  time  t  =  0.1  s  and  the  second  one  between  1000  and  800 
W/m^  at  time  t  =  0.2  s.  The  Figure  8  shows  that  the  PV  power  captured  varies  between  35.1  kW  and  71  kW 
and  then  returns  to  57.2  kW.These  values  correspond  is  shows  in  Figure  3  to  maximum  points  (M3,  Ml  and 
M2)  of  the  curves  associated  to  the  considered  radiation,  respectively.  The  Figure  8  also  shows  that  the 
voltage  of  the  PV  array  Vpv  varies  between  Vpv  =  648.8  V  and  Vpv  =  664  V  and  then  returns  to  661V,  which 
correspond  very  well  to  the  optimum  voltages.  Figure  9.  shows  the  injected  current  and  the  grid 
voltage  Vg^  This  Figure  9  clearly  shows  that  the  grid  current  is  sinusoidal  and  in  phase  with  the  grid 
voltage  Vg^,  proving  that  the  power  factor  unit  is  achieved.  The  alternating  currents  injected  to  the  grid  are 
illustrated  by  Figure  10. 
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(c) 


(d) 


Figure  7.  State  variables  and  their  estimates  in  transitional  regime  and  in  presence  of  radiation  step  changes 
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Figure  8.  MPPT  capability  of  the  controller  in  presence  of  radiation  step  changes 


Figure  9.  Unity  PF  achievement  in  presence  of  radiation  step  changes 
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Grid  citmcnls  ia.  ib  and  ic  [A] 


Figure  10.  A  zoomed  waveforms  of  the  alternating  eurrents  injeeted  to  the  grid 
in  presenee  of  radiation  step  ehanges 


4.2.  Temperature  variation  effect 

Figure  1 1  shows  the  performanees  of  the  eontroller  in  presenee  of  temperature  step  ehanges  while 
the  radiation  X  is  kept  eonstant  equal  to  lOOOW/m^.  The  temperature  step  ehange  is  performed  at  time  t=0.15s 
between  T=25°C  (298. 15K)  and  T=55°C  (328. 15K).  It  is  worth  noting  that  these  ehanges  are  very  abrupt 
whieh  is  not  realistie  in  praetieal  ease.  Nevertheless,  with  this  important  ehange  we  want  to  show  a  good 
robustness  of  the  proposed  eontroller  to  aehieving  the  MPPT  objeetive.  It  ean  be  seen  from  Figure  1 1  that  the 
system  traeks  the  new  operating  point  very  quiekly.  Indeed,  the  eaptured  PV  power  P  aehieves  the  values 
71.7kW  et  61.7kW  eorresponding  to  the  maximum  points  (M4  and  M5)  assoeiated  to  the  temperatures  25°C 
and  55°C,  respeetively  as  shown  in  Figure  4. 
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Figure  1 1 .  MPPT  and  DC  bus  voltage  behavior  in  presenee  of  temperature  step  ehanges 
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Figure  12  illustrates  the  grid  current  and  the  grid  voltage  Vga,  Figure  12  also  shows  that  the 
current  is  sinusoidal  and  in  phase  with  the  voltage  Vga  which  proves  the  unity  PF  requirement.  Finally, 
Figure  13  shows  the  zoomed  three  phase’s  grid  currents  and  voltages. 


Figure  12.  Unity  PF  behavior  in  presence  of 
temperature  changes 


Figure  13.  Zoomed  grid  currents  and  voltages  in 
presence  of  temperature  step  changes 
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5.  CONCLUSION 

We  have  addressed  the  problem  of  current  sensorless  control  of  the  photovoltaic  energy  conversion 
system.  The  main  feature  of  the  system  presented  is  it  does  not  require  an  intermediate  stage  of  DC-DC 
control,  because  the  maximum  power  is  set  by  the  inverter  itself  by  means  of  a  Perturb  &  Observe  algorithm. 

The  system  dynamics  has  been  described  by  the  nonlinear  state-space  model  (5a-b).  Based  on  a 
transformed  model  in  a- axis,  an  output  feedback  controller,  combining  an  observer  and  a  state  feedback 
control  laws  is  designed  and  analyzed  using  a  Fyapunov  approach.  The  controller  objectives  are  threefold:  i) 
ensuring  the  MPPT  in  the  side  of  PV  generator;  ii)  guaranteeing  a  power  factor  unit  in  the  side  of  the  grid,  iii) 
ensuring  the  global  asymptotic  stability  of  the  closed  loop  system.  Using  both  formal  analysis  and  simulation, 
it  has  been  proven  that  the  obtained  controller  meets  all  the  objectives.  It  is  worth  noting  that  all  objectives 
are  achieved  without  resorting  to  the  measurement  of  all  variables.  The  use  of  the  estimator  is  motivated  by 
the  global  cost  and  reliability  considerations.  Indeed,  the  fewer  the  number  of  sensors,  the  lower  the 
global  cost. 
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